Fast Condensation in a tunable Backgammon model 
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We present a Monte Carlo study of the Backgammon model, at zero temperature, in which a 
departure box is chosen at random with a probability proportional to (2ui — l)k + (1 — u) TV, where k 
is the number of particles in the departure box and TV is the total number of particles (equivalently, 
boxes) in the system. The parameter u) £ [0, 1] tunes the dynamics from being slow (u> — 1) to being 
fast (u) = 0). This parametrization tacitly assumes a two-box representation for the system at any 
instant of time and ui is formally related to the 'memory' parameter of a correlated binary sequence. 
For lu < 1/2, the system undergoes a fast condensation beyond a certain time that depends on ui 
and the system size TV. This condensation provides an interesting contrast to that studied with Zeta 
Urn model in that the probability that a box contains k particles evolves differently in the model 
discussed here. 

PACS numbers: 05.20.-y, 05.10.Ln, 02.50.Ey, 61.43.Bn 



Backgammon model [1-4] is a simple adaptation of the 
Ehrenfest's particles-in-boxes model, which demonstrates 
that the slow evolution of a system towards its ground 
state could be due to entropy barriers - sets of equal en- 
ergy states. Every box in this model is assigned a two- 
valued energy - namely, minus one if it is empty and zero 
if it is non-empty - so that the ground state of the system 
corresponds to having all the particles in one box. The 
system evolves from an initial zero energy state, corre- 
sponding to a " no empty box" configuration, towards its 
ground state by shifting a randomly chosen particle to a 
randomly chosen box subject to the standard Metropolis 
criterion. Since a particle from a more populated box 
is shifted out with higher probability, accumulation of 
particles in a box is a rare event, which is what is re- 
sponsible for the entropy barriers. A speed-up of the 
system dynamics can be achieved [5] by randomly choos- 
ing a non-empty box from which a particle can be shifted 
out. At zero temperature, the system ultimately reaches 
its ground state configuration, and density has no role to 
play in this condensation process. 

On the other hand, if an empty box is assigned zero 
energy whereas a non-empty box is assigned a negative 
energy proportional to the logarithm of the number of 
particles in it, then the system has been shown [6] to 
have a condensation transition at finite temperature for 
densities above a critical value. The process of conden- 
sation in this model, referred to as the Zeta Urn model, 
has been shown [7] to have a scaling behaviour and also 
an aging dynamics. 

In this context, it is of interest to know whether a 
parameter-driven condensation can be realized at all us- 
ing the Metropolis dynamics of the Backgammon model. 
Recognizing the formal analogy between this model and a 
binary sequence with long-range memory [8], we present 
Monte Carlo evidence to show that a tunable ('memory') 
parameter exists, below a critical value of which a con- 
densation transition can be realized at zero temperature. 
This study provides an interesting contrast to the the 
density-driven condensation of the Zeta Urn model. 



Two-box problem - The Metropolis dynamics of the 
Backgammon model, at zero temperature, has been rec- 
ognized [9,10] as equivalent to that of a two-box system 
consisting of a certain number of particles in one box 
and the remaining particles in the other box. If, at any 
point of time, the box containing less number of parti- 
cles is labelled zero and the other box labelled one, then 
the particles themselves can be identified by the labels 
assigned to the boxes in which they are found. The con- 
figurational state of this TV particles in two-box system 
can therefore be represented by a sequence of zeros and 
ones, {ai \ at =0,1;*= 1,2,..., TV}. This TV-tuple corre- 
sponds to a vertex of an TV-dimensional unit hypercube. 
A random walk on the vertices of this hypercube, in the 
presence of traps at (0, 0, ... , 0) and (1, 1, . . . , 1), then 
clearly represents the stochastic evolution of the two-box 
system. 

Thus, the Backgammon model may be viewed as a 
stochastic reduction or coarse-graining of the number 
of labels available for a particle to just one label ul- 
timately - i.e., it describes a stochastic procedure for 
coarse- graining a many-alphabets sequence into a single- 
alphabet sequence [11], after sampling enormous number 
of binary sequences en route. 

Let there be a total of TV particles distributed in two 
boxes such that no of them are found in a box labelled '0' 
and the remaining (n\ = TV — uq) are found in the other 
box labelled '1'. We assume that the box containing less 
number of particles is always labelled '0'. The probabil- 
ity, poj of choosing a particle from the box labelled '0' is 
given by (a la Ritort) 



which is also the probability of choosing the box labelled 
'0'. By adding and subtracting TV/ 2 in the numerator, 
we are led to a more general definition, 
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where n is a tunable parameter [8] . What this parameter 
stands for becomes clear if we rewrite the above definition 
in the form, 

p = w ^ + (l-u)^; W =(l + M )/26[0,1] (3) 

Now, ui may be interpreted as the weight associated with 
box zero. Since to lies in the unit interval, the parameter 
/i ranges from —1 to +1. The case /j, — 1 corresponds 
to choosing a box with a probability proportional to the 
number of particles in it, whereas the case ji = — 1 corre- 
sponds to choosing a box with a probability proportional 
to the number of particles in the other box. 

While the choice fi = 1 is known to slow down the sys- 



tem, the choice [i — — 1 on the contrary is expected to 
speed up the system; the parameter /i is therefore helpful 
in tuning the evolution of the system from being fast to 
being logarithmically slow. In particular, the mean relax- 
ation time for the system to reach the ground state from 
an arbitrary initial state should depend on the parameter 
/i. This parametric dependence as well as the nature of 
the transition from slow dynamics to fast dynamics can 
be studied by computing the Mean First Passage Time 
(MFPT) [10] for the general case of p given by Eq.(2). 

Using the same notation as in reference [10] , the MFPT 
Fjj—i(2L,fi) for a system, consisting of even number of 
particles (N = 2L) [12], to go from state j to state j — 1 
is given by 



L 

F j ,j-i(2L,n)=u j (2L, f i)+ u > 

n=j+l 



where 

«fc(2L,M) = \ ( 1 + ^ 2L 2L 2k] ) (6) 

We can easily compute them and show that they are 
monotonically increasing functions of L — j for all non- 
negative values of /i; i.e.,Fj j_i (2L, /i) > i^ +l j (2L, //) 
for all values of fi € [0, 1]. 

In fact, guided by the exact result F 10 (2L,/j, = 1) = 
2 2L ~ 1 — 1 for /i = 1, we have plotted in Fig.l the com- 
puted values of ln(Fi^(2L, fi) as a function of L for var- 
ious values of [i in the range < [i < 1. It is clear 
that Fi j0 (2L, /i) ~ a([i) L where the values of a(/z) for the 
/j, values considered are mentioned in the figure caption. 
Moreover, for [i = 0, we have the exact linear behaviour 
Fi, (2£,/i = 0) oc L. This implies that F li0 (2L,n) 
predominantly contributes to the system dynamics for 

Me [0,1]. 

However, for negative values of fi, the MFPTs do not 
have this monotonic behaviour. In Fig. 2, we have plotted 
Fj j_i(2L, n) as a function of j £ [2, 64], with L = 64, for 
various values of [i. Deviation from monotonicity may 
be observed even for a very small negative value of [i(= 
—0.005). Starting from the state j = L, the system slows 
down until it reaches the state, j c (A t ) < 64, after which 
its evolution is progressively fast. 

It is interesting to note that Fi i o(2L,/z) has the func- 
tional form e x l h for ji = —1.0 (Fig. 3). Since the prob- 
ability of choosing a box is proportional to the number 
of particles in the other box for the case fi = — 1.0, we 
expect that Fi.o(2L — > oo,/i = —1.0) — > 1, which is what 
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FIG. 1: Semi-log plot of the MFPT Fi,o(2L, /i) as a function 
of the system size L for /i = 1.0,0.8,0.6,0.4 and 0.2 (from 
top to bottom). The computed data suggest a functional 
form, Fifi(2L,fi) ~ a,(fj,) L , where a(fi) ~ 4,2.5,1.9,1.5 and 
1.2 respectively for the values of /x considered. 



we observe also in Fig. 3. For any other negative value of 
/i, Fi i o(2L, (j,) does not have a simple functional form. 

The above MFPT calculations clearly indicate that 
there is a speed up of the system's evolution for nega- 
tive values of fj,. In order to understand how the more 
populated boxes gain at the expense of the less populated 
boxes, we have carried out a Monte Carlo study of this 
problem. 

Tunable Backgammon model, Monte Carlo study - The 
configurational state of a system of N particles in N 
boxes at any given time is specified by the TV-tuple, 
{ni,n2, ■ ■ ■ ,tin}-, where is the number of particles 
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FIG. 2: A plot of Fj,j-i(2L,n) as a function of j € [2,64] 
for a system of size L = 64. The MFPT has been computed 
for fi = 0.005, 0.0, -0.005, -0.1, -0.3 and -1.0 (from top to 
bottom). 
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FIG. 4: Average number of particles per non-empty box, 
A(t), as a function of time. Data correspond to u = 
0.0, 0.25, 0.5, 0.75, 0.9 and 1.0, from top to bottom, for a sys- 
tem size N — 1024. Here, the actual number of Monte Carlo 
moves is equal to Nr. 
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FIG. 3: The MFPT Fi,o(2L,/x) computed and plotted as a 
function of the system size L for fi — —1.0. The slope of the 
straight line is ~ —1.03, suggesting thereby the functional 
form, Fi, (2L,/Lt = -1.0) ~ e 1/L . 



found in the k box. The evolution of the system is 
based on the dynamical rule that a non-empty box be 
picked up at random and a particle from it be trans- 
ferred into one of the remaining non-empty boxes. Thus, 
at every instant of time, there is a departure-box and a 
destination-box, which may be chosen as follows. 

We first assign to each of the non-empty boxes the 
normalized departure probability, dk = Pk/ ^ZuPki where 
Pk is defined in Eq.(3). We then use these departure 
probabilities to choose a departure-box at random. On 
the other hand, a destination-box is just one of the non- 
empty boxes chosen at random with equal probablity. 

Starting from the initial configuration, {np. = 1 | k = 
1,2,..., N}, we obtain the average number of particles 
per non-empty box, A(r), as a function of the Monte 
Carlo time r expressed in units of the system size N. In 
Fig. 4, we present the data obtained for a system of size 



N = 1024 corresponding to various values of ui. There 
is the expected [3] linear growth of A(t) with time for 
u> = 1/2. Linear growth is seen also for smaller values of 
uo < 1/2 upto a certain time t q (lu, N) beyond which the 
growth becomes superlinear (to(uj, N) — > oo as u> — > 0.5). 
Evidence for the existence of tq(uj, N) and the subsequent 
rapid growth is presented in Fig. 5. for the case ui = 0. In 
fact, A(t) saturates fast to the value N beyond t (uj = 

0, N). Quick estimates of tq(oj — 0, N), presented in the 
inset as a double logarithmic plot, suggest a power-law 
dependence t (lu = 0, N) ~ N 3 / 2 . 

On the other hand, for 1/2 < u> < 1, the growth of A(r) 
is linear upto a certain time tq beyond which it becomes 
sublincar (to — > as u> — » 1). It is a well established 
result [3] that A(r) oc ln{r) for the case u> = 1. In or- 
der to check whether the asymptotic sublinear growth 
seen for the other values of w(> 1/2) is also of this 
type, we present A(r) as a function of ln{r) in Fig. 6. 
The simulation data clearly show that the asymptotic 
growth for l/2<w<lis also described by the relation, 
A(t) oc ln{r). 

This may be anticipated from the departure proba- 
bility, dk, assigned to a non-empty box containing fc(= 

1, 2, . . . , N) particles: 

_ Pk = {2uj - l)k + {I - u)N 
k ~Z k Pk (2w-l) + (l-w)JV(l-/ ) U 

where /o is the fraction of empty boxes. In the initial 
stages of the system evolution (r < t (oj,N), for some 
To that depends on u and N), the number of particles 
in any given box will be much smaller than the system 
size (i.e., k << N). So, whatever be the value of ui, the 
dominating second term leads to 

dk - z —7- (8) 

1 - Jo 
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FIG. 5: Finite size effect on the average number of particles 
per non-empty box, A(t), for the case u = 0. There is a 
time, to(lj = 0,N), beyond which we see deviation from the 
linear growth to a rapid saturation at A(r) = N. Inset: A 
double logarithmic plot of the approximate estimates of to(N) 
for N = 128, 256, 512 and 1024. It suggests a power-law, 
to (AT) ~ N 3/2 . 

which is actually the uniform departure probability as- 
signed to the non-empty boxes for the case u> = 1/2. 

As the system evolves beyond tq(ui,N), there will be 
more boxes containing larger number of particles (i.e., 
k ~ N) and therefore, deviations from the behavior dic- 
tated by Eq.(8) are expected. That is to say, the first 
term in Eq.(7) progressively decides the system dynam- 
ics. For u> > 1/2, there will be a few highly populated 
boxes and the positive first term will ensure that the 
system evolution is similar to what we see for the case 
u> = 1. On the other hand, for u> < 1/2, the negative 
first term will ensure that the less populated boxes are 
emptied fast. 

Summary - We have shown that the standard Metropo- 
lis dynamics of the Backgammon model leads to a con- 
densation process if, at every instant of time, a depar- 
ture box is chosen at random with a probability given 
by Eq.(7) for the parameter u> < 1/2. While assigning 
this parametrized form of departure probability to a non- 
empty box, we tacitly assume a two-box representation 
for the system. The cases u> = 1/2, 1 correspond to Model 
B[3,5] and Model C[l-3] respectively. For all the other 
values of u, the system evolves like Model B till a cer- 
tain time tq(oj,N) beyond which it behaves like Model 
C for u> > 1/2 and, on the other hand, undergoes a fast 
condensation for uj < 1/2. That this condensation pro- 
vides an interesting contrast to that studied with Zeta 
Urn model becomes clear if we observe the temporal evo- 
lution of the probability that a box contains k particles 
(Fig. 7, for the case u = 0). Ultimately, the distribution 
is made up of two delta functions, one at k = and the 
other at k — N . 

I thank K. P. N. Murthy for helpful discussions. 
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